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Abstract. We present an implementation of Redelemeier's algorithm for the enumeration 
of lattice animals in high dimensional lattices. The implementation is lean and fast enough 
to allow us to extend the existing tables of animal counts, perimeter polynomials and series 

e expansion coefficients in rf-dimensional hypercubic lattices for 3 < < 10. From the data we 
, compute formulas for perimeter polynomials for lattice animals of size n < 11 in arbitrary 

dimension d. When amended by combinatorial arguments, the new data suffices to yield 
explicit formulas for the number of lattice animals of size n < 14 and arbitrary d. We also 
use the enumeration data to compute numerical estimates for growth rates and exponents in 
1^ ' high dimensions that agree very well with Monte Carlo simulations and recent predictions 

from field theory. 



rj^ • PACS numbers: 64.60.ah, 64.60.an, 02.10.Ox, 05.10.-a 



1. Introduction 



A polyomino of size n is an edge-connected set of n squares on the square lattice, a polycube 
of size n is a face-connected set of n cubes in the cubic lattice. Polyominoes and polycubes 
are a classical topic in recreational mathematics and combinatorics 1 1 J. In statistical physics 

■ and percolation theory, polyominoes and polycubes are called lattice animals Lattice 
\ animals are not restricted to dimension 2 or 3: a li-dimensional lattice animal of size n is a set 

of n face-connected hypercubes on If' . 
\^ \ In this contribution we will address the problem of counting the number affixed animals 

■ of size n in dimension d > 2. Fixed animals are considered distinct if they have different 
shapes or orientations. Free animals, on the other hand, are distinguished only by shape, not 
by orientation. Figure [T] shows all fixed polyominoes of size 3 and all free polycubes of size 
4. 

ITJi • We denote the number of cZ-dimensional fixed animals of size n by Aj(n). There is no 

[ formula for Aii{n), but we know that A^(n) grows exponentially with n. Using subadditivity 

■ and concatenation arguments IH, one can show that there are constants 1 < < 0° such that 

lim ^:/AAnj = Xd. (1) 

The constant A2 is known as Klarner's constant. 

A slightly stronger result due to Madras [5] asserts that 

lim . , , = Ad . (2) 
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Figure 1. All fixed polyominoes of size « = 3 (top) and all free polycubes of size n = A 
(bottom) and their perimeters. 



Intuitively, the growth rate should grow with the coordination number 2d of the lattice. In 
fact, in ||6l it is shown that 

Xci = 2de — o{d) , (3) 

and in the same paper it is conjectured that — {2d — 3)e + 0(1 /d). For finite d, however, 
we know only lower and upper bounds for Xj. Numerical estimates for Xj can be derived 
from extrapolating A J (n + l)/A^/(n), which is one motivation to compute A^/(n) for n as large 
as possible. We will try our hands at that in Section|6] 

In percolation theory one is interested in counting lattice animals of a given size n 
according to their perimeter t, i.e., to the number of adjacent cells that are empty (see 
Figure[T]i. If each cell of the lattice is occupied independently with probability p, the average 
number of clusters of size n per lattice site reads 

LgiVli-p)', (4) 

t 

where g'^^j denotes the number of fixed tZ-dimensional lattice animals of size n and perimeter 
t. The g's define the perimeter polynomials 

Pd{n,q) ^Y.gn tq' ■ 

t 

We can easily compute Aj(n) from the perimeter polynomial Pd{n,q) through 

A,in)=P,in,l)=Y,gi'l 

t 

In fact we can compute A j(« + 1) from the perimeter polynomials up to size n. 



(5) 



(6) 



(-1)"- 



(7) 



This equation follows from the observation that below the percolation threshold pc, each 
occupied lattice site belongs to some finite size lattice animal. 



n=\ 



(8) 
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for p < pc- The right hand side is a power series in p, and Equation (Q follows from the fact 
that the coefficient of p"^^ must be zero. 

As we will see in the next section, the algorithm for counting lattice animals keeps track 
of the perimeter anyway. Hence it is reasonable to use the algorithm to compute the perimeter 
polynomials and to apply (|7]i to get an extra value of A^/. 



2. The Algorithm 

The classical algorithm for counting lattice animals is due to Redelmeier [7|. Originally 
developped for the square lattice, Redelmeier's algorithm was later shown to work on 
arbitrary lattices and in higher dimensions ||8] and to be efficiently parallelizable ||9l. For two 
dimensional lattices there is a much faster counting method based on transfer matrices fTOl, 
but for d > 3 Redelmeier's algorithm is still the most efficient known way to count lattice 
animals. 

The algorithm works by recursively generating all lattice animals up to a given size nmax- 
Given an animal of size n, the algorithm generates animals of size n + 1 by adding a new 
cell in the perimeter of the given animal. The lattice sites that are available for extending 
the current animal are stored in a set U called the untried set. To avoid generating the same 
fixed animal more than once, lattice sites that have previously been added to the untried set 
are marked on the lattice. 
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Figure 2. Part of the square lattice that can be reached from animals up to size 4 (left). The 
number of lattice animals of size n <A equals the number of subgraphs in the neighborhood 
graph (right) that contain vertex 1 . 



In order to break the translational symmetry we demand that the initial site, which is 
contained in all animals, is an extremal site with respect to the lexicographic order of lattice 
coordinates. Figure |2] illustrates how this can be achieved in the square lattice. We simply 
block all lattice sites from further consideration that are in a row below the initial site (marked 
with a cross) or in the same row and to the left of the initial site. The generalization Xo d>2 
is straightforward. These blocked sites are never added to the untried set, but they need to be 
taken into account when we compute the perimeter f . 

We start with all lattice sites being marked "free" or "blocked," except for the initial site, 
which is marked "counted." Furthermore n—\,t — \ and the initial site being the only element 
of the untried set U . Redelemeier's algorithm works by invoking the following routine with 
this initial settings: 

• Iterate until U is empty: 
(1) Remove a site s from U. 
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(2) F :— set of "free" neighbors, B :~ set of "blocked" neighbors of s. N := \U\ + \B\. 

(3) Count new cluster: increase g„j+N-i by one. 

(4) If n < nynax- 

(a) Mark all sites in F and B as "counted." 

(b) Call this routine recursively with U'^UUF,n' = n+l and t' = t+N-l. 

(c) Relabel sites in F as "free" and sites in B as "blocked." 



Since the algorithm generates each lattice animal explicitely, its running time scales like 
Aii{n). This exponential complexity implies hard Umits for the accessible animal sizes. All 
we can do is to keep the prefactor in the time complexity function small, i.e., to implement 
each step of Redelemeier's algorithm as efficiently as possible by using an appropriate data 
structure for the untried set and for the lattice, see [8 1. 

Another crucial element of tuning Redelmeier's algorithm is the computation of the 
neighborhood of a lattice cell. In a recent paper [ 11] , Aleksandrowicz and Barequet observed 
that Redelmeier's algorithm can be interpreted as the counting of subgraphs in a graph that 
represents the neighborhood relation of the lattice. Figure [2] illustrates this for the square 
lattice. 

For any lattice, the neighborhood graph can be precomputed and be represented as an 
adjacency list which is then fed to the actual subgraph counting algorithm. That way the 
computation of the neighbors of a lattice cell is taken out of the counting loop, and the 
prefactor in the exponential scaling is reduced. 

The size of the neighborhood graph or, equivalently, the number of lattice points required 
to host lattice animals of size n scales like 0(«'^). Aleksandrowicz and Barequet ifTTl 
claimed that this exponential growth of memory with d represents a serious bottleneck for 
Redelmeier's algorithm in high dimensions. In a subsequent paper 1 12], they therefore present 
a variation of the algorithm that avoids the storage of the full graph by computing the relevant 
parts of the graph on demand. This cuts down the space complexity to a low order polynomial 
in d, but it forfeits the gain in speed that can be obtained by precomputing the complete 
neighborhood graph. 

We claim that in practice the space complexity of Redelmeier's algorithm is no 
bottleneck. The reason is that the prefactor in the &{n'^) scaling can be made small enough to 
hold the complete graph in memory for all values of n and d for which Af/(n) is computable 
in reasonable time. 

The key observation is that in [TT1, Aleksandrowicz and Barequet used a hyper-cube of 
side length 2n of the lattice to host the lattice animals, whereas a hyper-sphere of radius n 
suffices. This is a significant difference, as can be seen from the analogous situation in R''. 
Here the volume of a cube is much larger than the volume of the inscribed Euclidean sphere. 



Let B{d,n) denote the number of lattice sites that are n steps or less away from the origin. 
This "volume of the crystal ball" can be computed recursively via 

B(l,n) =2n + l 



• Return. 




(9) 



In the lattice I/', the number of lattice points in a cube of side length 2n + 1 is 

C{d,n) = {2n + lY'. 



(10) 



B{d,n)^B{d~l,n)+2Y^B{d-l,k) 



(11) 



k=0 
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Figure 3. Volume needed to cage animals of size n on the hypercubic lattice of dimension 
d = 9. Cubic {□) versus spherical cages (o)- 



which reflects the fact that the cry stall ball in dimension d can be decomposed into {d — 1)- 
dimensional slices whose diameter decreases with increasing distance from the central slice. 
The recursion (flU tells us that B{d,n) is a polynomial in n of degree d which can easily be 
computed, see A001845 to A001848 on oeis.org for the polynomials for d — 3, ... ,6. Note 
that B{d,n) can also be computed through the generating function lfl3l 

^^^-iBid,ny. (12) 

V ' 77=0 

The number of lattice sites needed to cage lattice animals of size n is very close to 
B{d,n — l)/2 for spherical cages and C{d,n — l)/2 for cubical cages. Figure |3] shows both 
numbers for the case d = 9. As you can see, d = 9 and n = 12 requires memory on the 
Terabyte scale if one uses cubic cages, but only a few Megabytes for spherical cages. 



3. Performance 



Our implementation of the Redelmeier algorithm consists of two programs. The first program 
computes the neighborhood graph of a specified lattice and writes this graph as an adjacency 
list into a file. The second program reads this file and computes the corresponding perimeter 
polynomials. The programs are written in C++ and can be downloaded from the project 
webpage ifTSl . 

When run on a laptop with an Intel® Core'"'^ 2 Duo CPU at 2 GHz, the program 
enumerates perimeter polynomials at a rate of roughly 2-10^ animals per second. This means 
that generating and counting one lattice animal and measuring its perimeter takes about 100 
clock cycles, which is reasonable for a program compiled from C++. 

At this rate, our laptop needs 35 days to enumerate the perimeter polynomials for d — 9 
and n < 11 (see Table |2]l. Computing the next perimeter polynomial (n = 12) would take 
more than three years. Note that according to Figure |3] the neighborhood graph for d — 9 and 
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Table 1. Range of perimeter polynomials and animal numbers in dimensions d > 3 that have 
been found by exhaustive enumerations. The old perimeter polynomials are from ISlandlTil. 
the old values of Arf(n) are from II II and 1121 . 
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Table 2. Number of lattice animals in the hypercubic lattice for d = 6 . . . 9 obtained by direct 
enumeration. New results in boldface, the numbers of smaller animals are from dUEll and 
references therein. Note that in [16], A%(n) and Ag(n) for n <9 were computed rather than 
enumerated by the same method that we will use in Section |4]to extend this table to « < 14 
and all values of d. 



n — 12 easily fits into the memory of a run-of-the-mill laptop. These numbers illustrate that 
for for all practical purposes, the bottleneck of Redelemeier's algorithm is time, not memory. 

Using a parallel implementation |9| that we ran on a Linux cluster with 128 Intel® 
Xeon® 3.2 GHz CPUs, or for the most demanding computations, on a SciCortex SC5832 
with 972 MIPS64 6-core nodes, we could extend the table of known perimeter polynomials 
and animal counts considerably, see Table [T] The new values for li < 5 are 

A3(19)= 651459315795897, 
A4(16)= 692095652493483, 
A5(14)= 227093585071305, 
A5(15) =3689707621144614. 

The numbers for 6 < < 9 are given in Table |2] the corresponding perimeter polynomials 
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can be found on the project webpage |fT5l . Before we evaluate the resuhs, we will discuss a 
combinatorial argument that allows us to extend the enumeration data considerably. 

Note that the most demanding computation in this paper was the enumeration of the 
perimeter polynomial for n — \A m d = 6. On a single core of a MIPS64, this enumeration 
would have taken 77 CPU years, on our Laptop from above it would still have taken about 7 
CPU years. In practice we used a parallel implementation that ran on many cores (and several 
different machines) such that no computation took longer than two weeks wall clock time. 

4. Proper Animals 

A lattice animal of size n can't span more than n — 1 dimensions. This simple observation 
allows us to derive explicit formulas for Aj{n) for fixed n. Obviously A^iil) — 1 and 
Aii{2) = d. A lattice animal of size n = 3 is either a one-dimensional "stick" with d possible 
orientations or "L-shaped" and spanning 2 out of d dimensions. Within these 2 dimensions 
there are 4 possible orientations for the L-shaped animal (see Figure [T]!, hence 



For n = 4, we have again the "stick" that lives in one dimension, 17 animals that span 2 
dimensions and 32 animals that span 3 dimensions: 



where DX(«,/) denotes the number of fixed proper animals of size n in dimension /. An 
animal is called proper in dimension d if it spans all d dimensions. Equation ( flj] ) is due to 
Lunnon ITTll . If we know A^(n) for a given n and d < t/max, we can use (fT3T l to compute 
DX{n,d) for the same value of n and all d < t/max, and vice versa. 

Since DX(n,/) — for / > «, Lunnon's equation tells us that A^/(n) is a polynomial 
of degree n — 1 in d, and since Ao(n) ~ for « > 1, it suffices to know the values 
Ai (n),A2(n), . . . ,A„_i (n) to compute the polynomial A^(n). From our enumeration data 
(Table[T]|, we can compute these polynomials up to Aj/(11), see Table|3] 

In order to compute A^/ ( 1 2 ) , we need to know A n ( 1 2 ) or equivalently, DX (12,11). The 
latter can actually be computed with pencil and paper That's because an animal of size 12 in 
1 1 dimensions has to span a new dimension with each of its cells to be proper In particular, 
its cells can't form loops. Hence computing DX(12, 11) is an exercise in counting trees. This 
is true forDX(«,« — 1) in general, so let's compute this function. 

The adjacency graph of a lattice animal of size n is an edge labeled graph with n vertices, 
in which each vertex represents a cell of the animal and two vertices are connected if the 
corresponding cells are neighbors in the animal. Every edge of the adjacency graph is labeled 
with the dimension along which the two cells touch each other 

In the case DX(n,« — 1), every pair of adjacent cells must span a new dimension. 
Therefore the corresponding adjacency graph contains exactly n — 1 edges, i.e. it is a tree, 
and each edge has a unique label. There are two directions for each dimensions that we 
represent by the orientation of the edge in the tree. Hence DX(n,n — 1) equals the number of 
directed, edge-labeled trees of size n, where in our context "directed" means that each edge 
has an arbitrary orientation in addition to its label. 





In general we can write 




(13) 
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Ad{2) =d 
Aj(3) ^2cP--d 

A,(6) = 216^4 + 23 W^ + ^d 

9604 . ,„^„,5 20651 . 14927 , 120107 . 827^ 
M^) = -^d^ - 1078^5 + — - — + — - —'^ 

, 262144 ,7 26624 132320 65491 ,4 1615991 113788 52589 , 

Ad(S) = d' <:rH d^ -\ d^ d^ -\ d 

' 315 5 9 3 90 15 42 

WON 118098 o ^,^,,,7 447903 . 511082 c 23014949 ,4 1522261 , 
^'^'^ =^/-26244rf^ + ^rf^ _,5 + ___,4 _,3 

38839021 2 30089 

, , 8000000 Q 2720000 = 14272000 7 11092360 . 239850598 

A^(IO) = d' d^ H d' rf" H 

''^ ' 567 21 27 9 135 

14606026 ,4 1067389643 , 42595493 2804704 , 

+ ^13^^ U^' + 

, 857435524 .n 67319318 ,0 2884481974 ,0 380707987 . 40341440233 . 

A^(ll) = d'" rf^H d^ rf'H d^ 

^ ^ 14175 105 945 45 2700 

1260803635 , 79118446751 4 19252021283 3 17126616179 , 7115086 

+ 5670 ' 252(^^ + 6300 ' 

^ 509607936 ,,, 15925248 .,) 607592448 ,n 1956324864 ,0 2930444704 ,7 

A^(12) = <i" rf"'H d^ rf**H d^ 

'^^ ' 1925 5 35 35 25 

2522387284 . 17894522696 , 1242881121 4 

rf" H <:r d^ 

15 105 10 

22272055467 , 4225468993 181356011 

+ d^ d^ H d 

350 210 66 

^ 551433967396 ^12 75047226332 166095324499 48436628461 9 

A^{13) = d d H d d 

' 467775 4725 1701 135 

49499551181119 o 1335959158369 ,7 248648897740349 . 25156285613453 ,< 

H d** rf'H d'' d- 

56700 900 136080 15120 

757565736903221 . 5607318230581 12648671104037 . 135165335 , 

H d^ d^ ^ d^ d 

680400 10800 83160 6 

, , 4628074479616 ^ 23612624896 „ 3309261190144 ,, 304034058496 .a 

^''^'^^= 868725 297 + 6075 ^5 ^ 

12648090831712 9 553376997376 g 758347226205724 7 2633038200122 g 

+ 2025 ^ 45 + 42525 ^ ^5 ^ 

98388569956577 , 2734657007119 4 11824147558382 , 560344373791 ,7 

H H d"* rf^ 

6075 270 2475 330 



(5 



97500388612 , 
+ 273 ' 



Table 3. Number of lattice animals of sizes 2, . . . , 14 in hypercubic lattices of dimension d. 
The polynomials for « < 11 have been obtained by direct enumeration and confirm those listed 
in |6|. Polynomials for 11 > 11 have been computed from enumeration data and known values 
of DX(n,n-Jt). 
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The number of vertex labeled trees of size n is given by n"^^, the famous fomula 
published by Cayley in 1889 1 18 1. The number of edge labeled trees seems to be much less 
known, at least it is proven afresh in recent papers like [19|. The following nice derivation 
is from 16 |. Start with a vertex labeled tree of size n and mark the vertex with label n as the 
root. Then shift every label smaller than n from its vertex to the incident edge towards the 
root. This gives an edge labeled tree with a single vertex marked (the root). Since the mark 
can be on any vertex, the number of edge labeled trees equals the number of vertex labeled 
trees divided by the number of vertices. According to Cayley's formula, this number is n"^^. 
And since each directed edge can have two directions, we get 

DX(n,«-l) = 2"-Si"-^ (14) 

This formula has been known in the statistical physics community for a long time ||20l . We 
used it to compute DX(12, 11) and then A^(12) (TableO. 

We can proceed further along this line. To compute Af/(13) we have to extend our 
enumeration data by A 1 2 ( 1 3 ),..., A g ( 1 3 ) or equivalently by DX ( 1 3 , 1 2 ) , . . . , DX ( 1 3 , 8 ) . What 
we need are formulas DX(n,n — k) for ^ > 1. 

For k > \, there is no longer a simple correspondence between edge labeled trees and 
proper animals. We need to take into account that there are edge labels with the same value, 
that the adjacency graph may contain loops, and that some labeled trees represent a self- 
overlapping and therefore illegal lattice animal. A careful consideration of these issues yields 

DX(«, n - 2) = 2""^ n"-^ (« - 2) (9 - 6« + 2n^) , (15) 

see ||6l for the derivation of (fT5t . 

For k > 2, the computation of DX(n,n — k) gets very complicated and is better left to a 
computer. In Appendix A| we show that 

DX(n,n - k) = 2"-2<^+i«"-2'^-igi(n) , (16) 

where gk{n) is a polynomial of degree 3k — 3. Hence we can compute g^ from 3A^ — 2 data 
points, like the values of DX(n,n — k) for n = k, . . . ,4k — 3. Our enumeration data suffices to 
compute g2 and §3 with this method, but not g4. 

However, there is a trick that allows us to compute gi^ from many fewer data points. The 
free energy 

/„ = ilogA,/(n) 
n 

has a well defined 1 /d expansion whose coefficients depend on n. If we assume that these 
coefficients are bounded in the limit n — > most of the coefficients in gi^ are fixed, and we 
only need to know k + I data points to fully determine gi^. See Appendix B for the details 
of this argument. In our case this enables us to compute gi^ up to = 7, see Table |4l and 
consequently A^/( 13) and A^/(14), see Table[3] 

We actually know all data to compute A(/(15) with the exception of the number A7(15). 
On our laptop, the enumeration of the missing number A7(15) would take about 80 years. On 
a parallel system with a few hundred CPUs this would still take several months, which is not 
out of reach. Computing the formula for Arf(16), however, is definitely beyond the power of 
our machinery. 

Before we turn our attention to the analysis of the enumeration data we note that 
Lunnon's equation (fT3] l has a corresponding equation for perimeter polynomials; 



^«.t-2(d-i)n- 
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g2(n) = («-2)(9-6;j + 2n2) 

Si (n) = ^^-^ (- 1560 + 1 122;j - 679«^ + 360n^ - 104n'* + \2rr') 
6 

g^(n) = (204960- 114302;j + 41527«^ - 17523;7 V7404«'* - 2930;j^ + 828;j'' - Ui.n' + %n^) 
6 

g^(n) = ^(-3731495040+ 1923269040n-535510740n2 + 150403080n^ -42322743^'* 
360 

+ 12397445n^ - 4062240)1*' + 1335320;j^ - 356232n'* + 62240^'' - 6000n'° + 240n' ' ) 

geW = ^^^(1785362705280-939451308048n + 248868418932«2 -56265094748«^ 
360 

+ 11984445891n'* -2448081038n^ +535284255n^ - 127651774;j^ + 33940138n^ 
-9580440^'* +2398912n'° -440688n" + 51856;;'^ - 3424n'^ +96^''*) 

SI (") = (- 15601775208155 1360 + 85 163968967728896;i - 225 1770497891913671^ 

45360 

+ 4585470174542376/7^ -851686123590540n'* + 146137469433102/7^ 
-2444108066052377*' +4148836864606/7^ - 747463726205/!*^ 
+ 149724735468/7'' - 33793043592/7"' + 8322494124/i" - 1946680944/!'^ 
+ 363148352/7'^ -47679184/7'"* +4019904/7'^ - 193536/7"' +4032/7'^) 

Table 4. Polynomials gi(/i) that appear in XyX(n.n — k) )16t . The polynomials g2 ^6 can 
be found as g^ o in Appendix 2 of 1211 . As far as we know, the polynomial gi has not been 
published before. See the Appendix of this paper for the method how to compute the . 



G\I denotes the number of proper t/-dimensional animals of size n and perimeter t. Since 

G^fl = for c/ > 77 — 1, we can write 

For a given value of n, ( fTSl l represents the perimeter polynomial for general dimension d. 
Our enumeration data allowed us to compute the G^^^} and hence the formulas ( fTSl l for /i < 1 1 
(see ifTSl for the data), extending the previously known formulas for /i < 7 1 14 1 . A computation 
of the next formula P^i{\2,q) requires the knowledge of the perimeter polynomials for d <\\ 
and n < 12. The enumeration for d = 11 and n = 12 alone would take ca. 38 years on our 
laptop. 

5. Mean Cluster Size 

From the perimeter polynomials we can compute moments of the cluster statistics like the 
mean cluster size 

S{p) = - L nVP./(n, 1-P)= Y.bAr)p' ■ (19) 

P n=\ r 

The coefficients of the series expansion read 

<..w = |»^I.S(,,;_„)(-ir'-" 

= (r+l)2j„(r+l)+£„2£gl">f ' )(-!)'+'-». (20) 

n=\ t V+f nj 
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12 


2 


30 


56 


90 


132 


3 


114 


320 


690 


1272 


4 


438 


1832 


5290 


12252 


5 


1542 


9944 


39210 


115332 


6 


5754 


55184 


293570 


1091472 


7 


19574 


290104 


2135370 


10159252 


8 


71958 


1596952 


15839690 


95435172 


9 


233574 


8237616 


113998170 


883192392 


10 


870666 


45100208 


840643170 


8258076192 


11 


2696274 


229502616 


6017266290 


76196541732 


12 


10375770 


1254330128 


44178511010 


710151162432 


13 


30198116 


6307973352 


315024296150 


6540805549192 


14 


122634404 


34574952952 


2307462163110 


60831844077672 


15 


327024444 


171364602736 






16 


1460721616 








17 


3347244554 








18 


17795165832 









r 


d^l 


d = % 


d^9 


1 


14 


16 


18 


2 


182 


240 


306 


3 


2114 


3264 


4770 


4 


24542 


44368 


74322 


5 


280238 


595632 


1146834 


6 


3210074 


8012384 


17720514 


7 


36394302 


107053424 


272530194 


8 


414610014 


1434259248 


4198328082 


9 


4685293438 


19125485024 


64487361906 


10 


53201681162 


255662267296 


991886672898 


11 


600207546946 


3405928921264 


1522687319670 


12 


6800785109594 


45466350310880 


233996383280898 


13 


76649757121000 







Table 5. Series coefficients of the mean cluster size S{p) = Y,f'rP^ in hypercubic lattices of 

r 

dimension d. New values in boldface, older values from 1 8] (d = 3) and 1141 (rf = 4, . . . , 7) and 
references therein. 



Since we can compute A^(r + 1) from the perimeter polynomials P^i{n,q) for n < r via 
we can also compute the series coefficients bd{r) from this set of perimeter polynomials. If 
we happen to know Ad{r + 2), we can get an extra coefficient through 

bd{r+l)={'- + 2)Ad{r + 2) 

.t..(-.-.)E.i?C,^„)'-""'- 

This formula can be derived by solving Q for Y.t8n,tt and plugging the result into (|20] l. 

We used ( |2TI ) to extend the table of coefficients (Table [S]) for c/ > 8, since here we know 
perimeter polynomials only up to « = 11, but the cluster numbers A^ up to n = 14. 
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Figure 4. Growth rate for d = 9. The symbols are ?uq{n) computed from Eqs. )23t . The 
coiTection exponent A = 0.58 and the hne are the result of a numerical fit to the three leftmost 
data points. 
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Figure 5. Exponent for d = 9. The symbols are dg{n) computed from Eqs. (23). The 
correction exponent A = 0.50 and the line are the result of a numerical fit to the three leftmost 
data points. 



6. Growth Rates and Exponents 

The cluster numbers Af/(n) are expected to grow asymptocically as 

b 

Ad{n) - CA,^ n"®'' (1 + -r- + coiTections), (22) 

where the exponents 0^/ and A are universal constants, i.e., their value depends on the 
dimension d, but not on the underlying lattice, while C and b are nonuniversal, lattice 
dependent quantities ll22l . The universahty facilitates the computation of 0^ for some values 
of d using field theoretic arguments. We know 03 = 3/2 12311241 . 04 = 11/6 ||25]| and 
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diuin 


MC 


CIlllIIl 




3 


2.12169 


2.1218588(25) 


1.489 


3/2 


4 


2.58750 


2.587858(6) 


1.796 


11/6 


5 


2.92254 


2.922318(6) 


2.113 


2.080(7) 


6 


3.17838 


3.178520(4) 


2.232 


2.261(12) 


7 


3.38403 


3.384080(5) 


2.357 


2.40(2) 


8 


3.55484 


3.554827(4) 


2.441 


5/2 


9 


3.70057 


3.700523(10) 


2.489 


5/2 



Table 6. Growth rates Xj and exponents 6^ obtained from extrapolating tlie enumeration data. 
Tlie columns marked MC contain values from large scale Monte Carlo simulations I27I28I . 



0j = 5/2 (the value for the Bethe lattice) for d > = the critical dimension for animal 
growth ||26l . 

The enumeration data for Af/(n) can be used to estimate both Xj and 0^/. For that we 
compute Xdiji) and 0d(n) as the solutions of the system 

\ogAd{n -k)^\ogC+{n- k) log A,/(«) - 0,/(n) log(n - k) (23) 

for A; = 0, 1,2. We need three equations to eUminate the constant logC. Growth rate Xd and 
exponent 0^/ are obtained by extrapolating the numbers Xd{n) and @d{n) to « ^ oo. From ( |22] | 
we expect that 

logA,,(«)-logA,/ + ^ (24) 

„A+1 

for large values of n. We used the data points Xd[n) for the three largest values of n to fit the 
parameters logA^/, b and A in (|24] |. A plot of \ogXd{n) versus n^^^^ (Figure |4]i then shows 
that the data points in fact scale like ( |24] |. The resulting estimates for logA^/ are listed in 
Table |6] They agree very well with the high precision values from large scale Monte Carlo 
simulations 027112811 . 

The same approach can be used to compute the exponent 0^. Here we expect 

@d{n)^@d + ^. (25) 

Again we used the data points 0£/(n) for the three largest values of n to fit the parameters 0^/, 
b and A. Figure |5] shows that @d{n) in fact scales like ( |25] l. The resulting estimates for 0^ 
(Table |6l) deviate from the Monte Carlo results and the exact values by no more than 3%. 

The estimates for Xd and 0^ based on the current known values of A^(n) are much more 
precise than previous extrapolations based on shorter sequences Arf(«), see II14II16L 

7. Conclusions 

We have seen that the memory requirements of Redelmeier's algorithm can be kept low by 
using hyperspherical regions of the lattice. Even in high dimensions, the limiting resource in 
Redelmeier's algorithm is time, not space. 

We used a lean and efficient implementation of Redelmeier's algorithm to compute new 
perimeter polynomials in hypercubic lattices of dimensions d < 10. We have used these new 
perimeter polynomials together with combinatorial arguments based on proper animals to 
compute new values of the cluster numbers Aj(n) and new formulas for Ad{n) for n < 14 and 
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9 



9 




Figure Al. A labeled spanning tree that contains a part O — > O < — > O < — O corresponds 
to a 4-loop in the lattice animals, a quadrilateral that lives in the ('- j plane (left). To count the 
number of spanning trees with such a 4-loop, the edges of the quadrilateral are removed and 
the vertices of the quadrilateral are considered the root vertices of disconnected trees (right). 
The number of the latter is given by jA. U with i = 



arbitrary d. We've also used the new data to compute formulas for the perimeter polynomials 
Pci{n;q) for n < 11 and arbitrary d. We haven't shown these formulas here, but you can 
download them from the project webpage [ 15 1. 

We've also used our data to compute the formula for DX(n,n — 7), the number of proper 
animals of size n in dimension n — 7, and new coefficients in the series expansion of the mean 
cluster size S{p). 

Based on the enumeration data, we've finally computed numerical values for the growth 
rates and the critical exponents Od that agree very well with the results of Monte Carlo 
simulations and field theoretical predictions. 

All in all we have explored the limits of computerized counting of lattice animals in 
dimensions d > 3. Any significant extension of the results presented here would require 
either a considerable amount of CPU time or an algorithmic breakthrough comparable to the 
transfer matrix methods for d = 2. 

Appendix A. Formulas for DX(n,n — k): Structure 

In the physics literature like ||2T|, equation (fT&t is usually assumed to be true just because it is 
supported by the available enumeration data. But as a matter of fact, one can actually motivate 
(fTSI ) using the type of arguments that were used in |6 | to prove the formula for DX(n,n — 2). 
The idea is to show that the leading order of DX(«,« — k) is ^ 2"n"+*^^ whereas the lowest 
order contributions are ^ 2"n"^^^^'. This is exactly the range of terms in (fTST i if gk is a 
polynomial of degree 3(^— 1). 

Equation ( fTSI l is obviously correct for k — I (with gi — 1). For k > 1, we can still 
represent lattice animals by trees, namely the spanning trees of their adjacency graph. Each 
spanning tree is again an edge labeled tree, but this time there are only n — k labels for « — 1 
edges, i.e., k—l edges will carry a label that is also used elsewhere in the tree. Consider a tree 
whose n — 1 edges are labeled with numbers 1 , . . . , n — 1 . If we identify each of the high value 
labels n — k+l,n — k + 2,...,n— 1 with one of the low value labels 1, . . . ,n ~ k, we get the 
right set of labels. Since there are (n — k)^^^ ways to this, the number of edge labeled trees 
with n — k distinct labels scales like n'^^'n"^^ = «"+'^^^ in leading order With two directions 
for every edge we get 2"n"^''^'^ for the leading order in DX(«,n — k). 

For k > 1, a proper animal can contain loops. Figure IaTI shows the simplest case: a 
loop that arises because one dimension (/) is explored twice. The result is a quadrilateral that 
lives in the /- j plane of the lattice. In the spanning tree, this corresponds to a part that is 

labeled O — ^ O O O. Since graphs with loops have several spanning trees, our 
count of edge labeled trees overcounts the number of proper animals. We need to subtract 
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4-1 k-l i\ 

9 — <^ — *o o 

; j 

Figure A2. An animal that contributes to DX(«,« — k) uses at most k — \ dimensions more 
than once, and the longest loop arises when each of these k—\ dimensions is explored twice. 
Such a loop contains 2k vertices, as shown here. 



some contributions from loopy animals. The idea is to break up the part of the spanning tree 
that corresponds to the loop and to separately count the number of trees that are attached to 
the vertices on the loopy part. 

Consider an animal with a loop that contains I cells and one of its spanning trees. If we 
remove all the edges from the spanning tree that connect the vertices in the loop, the remaining 
graph is a forest, i.e., a collection of trees, where each tree is rooted in one of the £ vertices. 
The forest has n — i vertices with a total of n — £ edges, and it is independent from the way 
that the root vertices have been connected in the loop. 

Now the number of ordered sequences of £>l directed rooted trees with a total of n~£ 
edges and n — £ distinct edge labels is 

2"-^n"-'^-^£. (A.l) 

See ||6l for a proof of (lA.lb . The lowest order corrections come from those animals for which 
the number £ of cells in a loop is maximal. This is the case for a loop that joins all — 1 
non-unique edge labels, see Figure |A2] The number of vertices in these loops is £ = 2k, hence 
the lowest order corrections are ^ 2"^^^n"^^^^^ . 

If we want the exact number of spanning trees for loopy animals, we need to count the 
number of ways to reconnect the roots of the forest to form a single tree. But since this number 
does depend on k but not on n and we are interested only in the scaling with «, we don't need 
to enter this discussion here. The same is true for animals that contain several small loops 
instead of a single loop of maximal length. If we apply the separation trick to one of the 
shorter loops, we get a scaling of order larger than ^ 2"^^'^n"^^^^^ , and the resulting forest 
is then labeled with fewer labels than edges, which increases the order even further So the 
lowest order corrections from loops come in fact from single loops of maximum length, and 
these contributions are of order ^ 2"^^'^n"^^*^^', as claimed in (fTST l. 

Besides the non-uniquess of spanning trees for loopy graphs, there is another type of error 
that needs to be corrected: some edge labeled trees correspond to animals with overlapping 
cells, i.e., to illegal animals. For instance, if a spanning tree contains the subtree 

• ^0^# or ©^O-^*, 

the two ©'s represent the very same cell of the animal. But these "colliding" configurations 
can be interpreted as 3-loops, or more generally, as £ loops, and counted in the same way as 
the legal loops above. Again the lowest order contributions come from the longest "colliding" 
loops which are formed by A: — 1 labels assigned to two edges each and arranged like 

• -^o~^o---o-^o-hLo-^^o---o<^#. 

These longest collision loops contain 2k — 1 vertices of the tree (representing 2k — 2 cells of 
the animal). According to ( lA.lb . their number scales like 2"^2A;+i^n-2i:^ order above the 
lowest order of legal loops. 
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This concludes the motivation of ( fT6] l. Note that a proof of ( fT6] ) would require a thorough 
analysis to exclude contributions outside the range covered by ( fT6] l. 



Appendix B. Formulas for DX(n,n — k): Coefficients 

Having established the fact that DX(n,n — fc) is given by ( fT6l l we still have to determine 
the coefficients of the polynomials gk{n)- Since has degree 3(^—1), it seems that we 
need to know the 3k ~ 2 values DX{k, 0) , DX(A; + 1 , 1 ) , . . . , DX(4A: - 3, 3A: - 3) to compute the 
coefficients. In terms of our enumeration data this means knowledge of Arf(n) for n <4k~3 
and d <n — k. The data in Table |2] suffices to compute the coefficients of g/^ for k^2,3, but 
not for k>4. Nevertheless we can compute gi^ for A: < 7 by assuming that the "free energy" 

Hm - logA,;(n) 

has a well defined 1 /d series expansion. This approach has been used to compute g/^ for k <6 
from much less enumeration data in lISTI and ll29l . and we used it to compute gj from the new 
enumeration data. Since the method hasn't been described in detail elsewhere, we provide a 
description in this Appendix. 

Let's start with Lunnon's equation, which tells us that A(/(n) is a polynomial of degree 
n—1 ind with coefficients that depend on n. For c/ > « we have 



n-l 



A^(«)= £DX(n,^) 



k=l 



n-l 

= ^ aj{n)d' 

.7=1 



with 



where 



aj{n) 



1^DX{n,k) 



(B.l) 



(B.2) 



denotes the Stirling number of the first kind. In particular we get 

DX(n,n-l) 



fln-i(n) = 
a«-2(n) 
flH-3(n) = 



(«-l)! 
DX(n,n-2) 

(«-2)! 
DX(n,n-3) 



DX(«,«-1) 

DX(n,«-2) 
(«-2)! 



n-l 
n-2 
n-2 
n-3 



DX(n,n-l) 



(«-3)! 

and so on. From ( IB. Il l we get 

Arf(«) =fl„_i(ny""' 
and with ln(l +x) — x + x^/2 +x^/3 . . . this gives the 1/d series for the "free energy" 



n-l 
n-3 



•1-2 



1 



1 

-lnAj(«) 
n 



1 l^ 1 J , 1 1 f \ , ^ an-2{n) 1 

1 ln<i + -lnfl„_i(n) H - + i 

n J n n a„_i[n) d 



We assume that all coefficients in this series remain bounded in the limit n 
definitely true for the zeroth order term: 

lim -lna„_i(n) = 1 +ln2. 

n-Joo n 



°o. This is 
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For the first order coefficient we get 



1 fl„_2 n— lDX(n,n — 2) 1 
nfl„_i n DX(n,n— 1) n 



n-1 
n-2 



= ( 1 _l^ •^^(n) n-2 



n 



4n2 2 

This is only bounded if the g2{n) term balances the second term, i.e., if the n^ coefficient of 
g2 equals 2. Using also the fact that DX(2,0) = 0, we can write 

^2 (n) = [n — 2) (2n^ + bn + c) . 

To compute the remaining coefficents, we only need to know DX(3, 1) = 1 and 

DX(4,2)=A2(4)-2=17 

to get 

giin) = (n-2)(2«2-6« + 9). 

The postulation of bounded coefficients in the series ( IB.3l l has saved us from knowing the 
value DX(5, 3) to compute g2- How much does it help us to compute g^l 
The polynomial gi^ enters the series expansion (IB. 3) via the term 



1 DX(n,«-A:) («-l) 
nDX(n,«- 1) (n-fe)l 



= 2^-^^n 1 -^''gkin) {n-l){n-2)---{n-k+l) 



0(n*-l) 

in the coefficient of d^^'^^^^ The leading order of this term is n^^gi^{n). All terms of degree 
larger than k in the polynomial gi^ lead to unbounded contributions to the series coefficient 
that need to be counterbalanced by other terms. These balancing terms always exist, a 
fact that gives additional support for the claim of bounded coefficients. The coefficients 
of the terms of order larger than k in gi^ are therefore computable from the known terms 
gk-\{n),gk-2{n)T . . that also enter the same coefficient. Only the k+\ low order terms 
of gk are not fixed by the postulate of bounded coefficients and we need ^ + 1 data points 
DX(A:, 0) , DX(fe + 1 , 1 ) , . . . , DX(2/fc, k) to complete gk. 

Our enumeration data suffices to compute g^ (see Table |4]i. The computation of 
requires knowledge of DX(15,7) and DX(16,8), or in terms of A(/(n), 

DX(15,7) =A7(15) -572521427068702741 

and 

DX(16,8) =A8(16)+48366334433679758-56*A5(16) 
+ 28*A6(16)-8*A7(16). 
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